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1.  Introduction 


Parabolic  equation  (PE)  methods  have  become  very  popular  for  the  calculation  of  sound 
transmission  along  ground-to-ground  paths  in  the  atmosphere  [1].  These  methods  can  account  for 
refraction  of  the  propagating  waves  by  atmospheric  wind  and  sound  speed  gradients,  diffraction 
into  shadow  regions,  and  forward  scattering  of  sound  by  atmospheric  turbulence.  With  the 
introduction  of  the  Green’s  function  PE  by  Gilbert  and  Di  [2],  it  became  possible  to  perform 
sound  field  calculations  with  parabolic  equations  very  rapidly.  Nevertheless,  even  the  Green’s 
function  PE  is  too  slow  in  many  Army  applications,  where  calculations  are  required  at  many 
frequencies  in  fractions  of  a  second.  Among  these  applications  are  sensor  platforms  with 
algorithmic  adaption  to  the  propagation  environment,  tactical  decision  aids,  and  combat 
simulations. 

To  address  the  Army’s  need  for  faster  sound  transmission  calculations,  we  are  developing  an 
artificial  neural  network  model  to  predict  the  transmission  loss  (diminishment  in  sound  energy 
resulting  from  propagation  in  the  atmosphere).  If  successful,  the  neural  net  will  provide  an 
extremely  fast  but  still  accurate  determination  of  transmission  loss.  Although  neural  nets  have 
been  used  previously  in  acoustics  [3],  [4],  [5],  the  main  application  has  been  automatic  target 
recognition  (ATR).  Here,  we  are  attempting  something  completely  different,  namely,  the 
prediction  of  propagation  effects.  In  principle,  the  propagation  neural  net  could  be  used  to 
compensate  signatures  for  the  propagation  effects,  thereby  improving  robustness  of  ATR 
algorithms,  whether  based  on  neural  nets  or  other  methodologies. 

In  the  development  of  any  neural  net,  an  appropriate  set  of  input  (training)  and  output  variables 
must  be  determined.  The  output  variable  for  the  propagation  neural  net  is  simply  the  sound 
pressure  at  the  receiver  position.  The  input  variables  and  parameters  would  be  the  propagation 
geometry  (target  and  sensor  heights,  horizontal  separation,  sound  frequency),  atmospheric 
conditions  (wind  velocity  and  sound  speed  profiles),  and  acoustical  ground  properties. 

Successful  development  of  the  neural  net  requires  the  output  variable  to  have  a  fairly  smooth 
dependence  on  the  input  variables.  It  is  also  desirable  to  keep  the  input  variables  as  few  in 
number  as  possible.  For  this  purpose,  we  develop  in  this  report  a  nondimensional  version  of  the 
acoustic  propagation  model  that  describes  the  propagation  with  a  minimal  number  of  parameters. 
The  propagation  model  includes  the  wave  propagation  equation  (in  this  case,  the  PE)  as  well  as 
equations  for  the  boundary  conditions,  initial  condition,  and  atmospheric  fields. 


2.  Derivation  of  the  Nondimensional  Equations 


Among  the  features  of  this  model  are  input  parameters  that  specify  the  acoustic  properties  of 
porous  media  using  relaxation  theory  [6],  the  ability  to  handle  both  narrow  and  wide  angle  cases, 
attenuation  of  the  reflection  at  the  upper  boundary,  and  a  second  order  accurate  ground  boundary 
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condition.  The  model  was  written  in  MATLAB1  [7]  and  follows  the  procedures  that  are  outlined 
in  West  et  al.  [1],  which  use  a  finite  difference  numerical  technique  to  solve  the  parabolic  partial 
differential  equation. 

2.1  Parabolic  Equation 

We  take  as  our  starting  point  the  following  narrow  angle  parabolic  equation  for  sound 
propagation  in  a  moving,  inhomogeneous  atmosphere: 

dP/dx  =  i/(2k0)(d2/dy2+d2  /dz2  +k2s  - k2)p,  (1) 

in  which  p  =  P  ex p  (ikax)xm  is  the  sound  pressure,  x  is  the  horizontal  coordinate  in  the 
propagation  direction,  y  is  the  horizontal  coordinate  perpendicular  to  the  propagation,  -  is  the 
height,  ifceff  is  the  effective  wave  number,  and  ko  is  the  wave  number  at  ambient  conditions.2 
Furthermore,  in  which /is  the  frequency  and  ceff  =  c  +  ux  is  the  effective  sound 

speed,  c  being  the  actual  sound  speed  and  ux  the  component  of  the  wind  velocity  in  the  direction 
of  propagation.  This  equation  includes  waves  traveling  outward  from  the  source  in  a  roughly  30- 
degree  beam.  It  was  derived  by  Nghiem-Phu  and  Tappert  [8]  and  is  equivalent  to  Equation  (2.88) 
in  Ostashev  [9]  if  the  components  of  the  wind  perpendicular  to  the  propagation  are  neglected. 

The  PE  code  we  use  in  this  project  provides  a  “pseudo  three-dimensional”  solution  to  Equation 
(1),  as  described  by  West  et  al.  [1].  In  actuality,  the  field  is  determined  only  in  the  vertical  plane 
of  the  x-  and  z-  coordinates.  This  is  accomplished  by  neglecting  P/x  in  comparison  to  dP/dx  (a 
far-field  approximation)  and  dropping  the  differentiation  in y,  with  the  result 

dP / dx  =  i/ (2 k0 )(d2  /dz2  +  k2s  -  k2 )p.  (2) 

This  is  equivalent  to  Equation  (15)  in  West  et  al.  [1]  (although  with  /-it  replacing  k)  and 
represents  the  standard  narrow  angle  parabolic  equation.  Note  that  P  in  this  equation  does  not 
have  units  of  pressure  but  rather,  units  of  pressure  divided  by  the  square  root  of  length.  We  will 
refer  to  P  as  the  surrogate  sound  pressure.  In  practice,  we  may  set  £eff  =  2nflc&  +  ia  in  the 
preceding  equation,  in  which  a  is  an  attenuation  coefficient  that  represents  the  conversion  of 
sound  energy  to  heat  in  the  propagating  waves. 

Equation  (2)  contains  three  fundamental  physical  dimensions:  time,  length,  and  mass.  To 
develop  a  nondimensional  version  of  this  equation,  we  must  normalize  by  three  quantities  that 
involve  these  dimensions.  We  choose  here  the  ambient  sound  speed  (fo),  ambient  air  density 
(/d),  and  frequency  (/).  Indicating  normalized  variables  with  overbars,  we  have  the  following 
nondimensional  quantities: 

P  =  P{co  / /)' 2  /(/Vo ) 

x  =  xf  /  c0  (4) 

z  =  zf/c0  (5) 

1  MATLAB  is  a  registered  trademark  of  the  MathWorks. 

2 Ambient  in  this  report  refers  to  an  arbitrary  reference  value  that  is  characteristic  of  the  propagating  environment.  Typically, 
the  value  at  the  source  height  is  used. 
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*0  =  Kco  f  f  = 


(6) 


and 

^eff  =  /  f  (T) 

Substituting  these  quantities  into  Equation  (2)  results  in  the  following  nondimensional  version  of 
the  parabolic  equation: 

dP/dx  =  i /(4?r)(d2  /dz2  +  k2s  - 4k2 )p.  (8) 

This  equation  is  solved  in  a  manner  similar  to  the  procedure  used  for  the  dimensional  version  in 
West  et  al.  [1],  and  the  detailed  steps  can  be  obtained  from  this  reference.  First,  the  terms  are 
rearranged  so  that  the  equation  is  written  as 

8P/8x  =  ad2P  /  dz2  +  bmP  (9) 


and  the  coefficients  a  and  b„„  after  it  is  transformed  into  a  finite  difference  equation,  are  found  to 
be 


a  =  i /(4?rAz2) 


for  m  =  1,2, . . . ,  M 


(11) 


in  which  Mis  the  total  number  of  vertical  mesh  points  and  A z  is  the  nondimensional  difference 
between  two  vertical  mesh  points 

^  ~  Zm+ 1  —  Zm  (I2) 

The  complex  wave  number,  km ,  is  defined  by 

km  =2x/ceffm  +iam  (13) 

In  this  equation,  ceff  is  the  effective  sound  speed  in  layer  m  normalized  by  Co  and  am  is  the 
attenuation  normalized  by  wavelength  (cq/J). 

Once  P  has  been  determined,  the  actual  nondimensional  pressure  follows  from  the  equation 

p  =  Pe’2™/y[*.  (14) 

2.2  Ground  Boundary  Condition 

As  discussed  by  West  et  al.  [1],  the  ground  impedance  boundary  condition  for  the  parabolic 
equation  is 

dP/dx  =  -2mfp0P  /  Zc,  (15) 
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in  which  Zc  is  the  ground  impedance.  Many  models  are  available  for  the  impedance  of  porous 
ground  surfaces  such  as  are  encountered  outdoors.  A  particularly  simple  one  that  nevertheless 
gives  accurate  results  was  developed  by  Wilson  [10].  The  impedance  equation  is 


Poco 


Q 


1  + 


/-I 


V 


1 


(l-ia»r,f5X  (l  -i(orvf5 


-0.5 


(16) 


in  which 

rv  =  PS1  /(2m;Q)and  xe  =  Np[xv  (17) 

In  these  equations,  the  xv  is  a  relaxation  time  associated  with  viscous  diffusion,  xe  is  a  relaxation 
time  associated  with  thermal  diffusion,  y  is  the  ratio  of  specific  heats,  Npr  is  the  Prandtl  Number, 
q  is  the  tortuosity,  Q  is  the  porosity,  sp  is  the  pore  shape  factor,  and  cris  the  static  flow  resistivity. 
We  would  like  to  nondimensionalize  the  impedance  equation  in  a  form  using  fewer  parameters. 
To  do  so,  we  first  note  that 

i(ozv  =  ilcop0q2f  /(l os2pCl)  (jg) 

and  define  a  nondimensional  flow  resistivity,  a ,  which  depends  on  parameters  previously 
defined. 

ct  =  otJ/(/a>Q)  (19) 


Combining  Equations  (18)  and  (19)  gives 

id) xv  =  i7xq2  /(oQ.)2  and  icoxe  =  i7xNVJq 2  /(oQ.2 ) 


(20) 


Substituting  them  into  Equation  (16)  gives 


Zc  =  q\ 


1  + 


Y- 1 


(l -inNprq2  /  a) 


(l-irij2  /af5 


-0.5 


(21) 


in  which  Zc  =  Z/(/^,c0)  and  q  =  q/Cl.  Considering  that  y  and  Npr  are  determined  once  the 
temperature  is  specified,  the  only  parameters  that  are  required  to  calculate  Zc  are  cr  and  q  . 
Thus,  the  number  of  ground  parameters  has  been  reduced  from  four  (cr,  q,  Q,  and  sp)  to  two  (cr 
and  q). 

In  nondimensional  form,  the  boundary  condition.  Equation  (15),  is 

dP/8z  =  -2mP/Zc.  (22) 


Regarding  the  numerical  implementation  of  dP  /  8z  in  spatially  discretized  form,  the  first  order 
accurate  case  can  easily  be  obtained  (West  et  al.  [1],  Section  5).  For  the  nondimensional  case,  the 
first  order  accuracy  boundary  condition  is  found  to  be 
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(23) 


P0  =i>(l-2/nAz/Zc)_1 

and  the  respective  nondimensional  boundary  condition  for  second  order  accuracy  is 

Po  =  (4 Pi  -  p2  )(3  -  4/nAz  /  Zc  )■'  (24) 


2.3  Upper  Absorbing  Boundary  Condition 

In  order  to  avoid  spurious  numerical  reflections  of  the  sound  wave  at  the  upper  boundary  of  the 
computational  domain,  an  absorbing  layer  at  the  top  of  the  domain  is  often  added.  Salomons  [11] 
suggests  the  following  quadratic  form  for  the  attenuation  coefficient  in  this  layer: 


a  =  ju\ 


z-z. 


abs 


VZtop  ^abs  J 


z  >  z. 


abs 


(25) 


in  which  ju  is  an  adjustable  parameter  (dimensions  of  inverse  length),  za bS  is  the  height  where  the 
absorbing  layer  begins,  and  ztop  is  the  top  of  the  computational  domain.  The  absorbing  layer  is 
typically  30  to  100  wavelengths  deep.  Thus,  we  might  rewrite  Equation  (25)  in  nondimensional 
form  as 


a  =  /u 


V 

^abs 

V  J 


,  z>  z. 


abs 


(26) 


in  which  Nx  is  the  depth  of  the  layer  in  wavelengths,  a  =  aco/f,  and  //  =  /ucjf  Salomons  [11] 
suggests  values  for  the  parameter  //that,  while  frequency  dependent,  are  only  roughly 
proportional  to  f.  Here,  we  propose  setting  Ji  =  1 .  This  leads  to  values  of  //  that  are  fairly  close 
to  the  ones  in  Salomons’  Table  I,  particularly  near  100  Hz. 


2.4  Initial  Condition 

The  “initial  condition”  here  refers  to  the  sound  pressure  at  x  =  0,  from  where  the  solution  is 
marched  forward  in  distance,  as  opposed  to  the  solution  at  some  initial  start  time.  (Keep  in  mind 
that  the  source  is  assumed  to  be  harmonic,  so  that  there  is  no  start  time.)  Following  Salomons 
[11],  we  use  the  following  Gaussian  starter  function  for  the  initial  condition: 


P(0,z)=sjkl- [exp 


+  Rex  p 


in  which  zs  is  the  source  height,  S  is  the  source  strength  (in  units  of  pressure)  and 

Zc/{p0c 0)-l 


R  = 


Zc/{p0c  0)+l 


(27) 


(28) 
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(29) 


is  the  surface  reflection  coefficient.  In  nondimensional  form,  these  equations  become 
p(0,  z)  =  S  V2^r  jexp[-  2nl (z  -  zs  )2  J+  Rex p[-  In1  (z  +  zs  )2  JJ , 


in  which 


and 


Zc+1 


(30) 


(31) 


2.5  Atmospheric  Profiles 

By  the  refraction  of  sound  energy,  the  wind  and  temperature  profiles  near  the  ground  can  have  a 
very  significant  impact  on  the  pressure  field.  Therefore,  we  also  need  a  nondimensional  model 
for  these  profiles.  Let  us  begin  by  writing  the  effective  sound  speed  as 

ceS=c  +  ucos/3,  (32) 


in  which  c  is  the  actual  sound  speed,  u  is  the  horizontal  wind,  and/? is  the  angle  between  the 
propagation  direction  and  the  wind.  Neglecting  the  effect  of  humidity  (which  is  usually  small), 
the  sound  speed  in  air  is  given  by 


1  +  - 


r 


27; 


0  J 


(33) 


in  which  T  =  T0  +  T'  is  the  temperature  and  T  a  small  perturbation  about  the  ambient  value  T0. 
The  following  equations,  based  on  Monin-Obukhov  similarity  [12,  13],  are  known  to  describe 
the  near-ground  temperature  and  wind  profiles  well  over  reasonably  uniform  terrain. 


T{z)-T{zr)  = 


P,T, 

ln-^--  T/( 

f  _  \ 

+n 

^  /■ 

^  r 

V  L(>  ) 

V  Lo  )_ 

(34) 


and 


V 

(  V 

5“ 

1 11 
i 

- 

U1? 

w 

z0 

ky 

--0 

K^o  j 

V  ^ o  )_ 

(35) 


in  which  kv=  0.40  is  von  Karman’s  constant,  Pt  =  0.95  is  the  turbulent  Prandtl  number,  Yd  =  g!Cp 
=  0.0098  K/m  is  the  dry  adiabatic  lapse  rate  (accounting  for  the  decrease  of  temperature  with 
height  because  of  compression  in  the  air  column),  g  is  gravitational  acceleration,  z0  is  the  surface 
roughness  length,  u,  is  the  friction  velocity,  T,  =  -  <  WT  >s  /  u,  is  a  temperature  scale,  and 
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<  w'T'  >  v  is  the  covariance  of  vertical  velocity  and  temperature  at  the  surface  (kinematic  heat 
flux).  The  variable  zr  is  a  reference  height  at  which  the  temperature  is  measured;  if  we  explictly 
set  To=T ( zr •),  then  T  in  Equation  (33)  equals  T(z) -  T(zr).  The  T’s  are  universal  functions  of 
the  dimensionless  ratio  C,  =  z!L0,  in  which  La  =  u2tT0  /  kvg  <  w'T'  >(  is  called  the  Obukhov  length. 

The  following  equations  for  H'/,  and  T,,,  are  recommended  based  on  Hogstrom  [14]  and  Wilson 
[15]: 


2  In 


1  + 


/  2 


,£<0 

C>o’ 


(36) 


in  which  the  a’s  and  b’ s  are  constants  with  the  values  a*  =  7.9,  a„,  =  3.6,  by,  -  8.4,  and  bm  =  5.3. 
The  case  £<0  corresponds  to  buoyantly  unstable  conditions,  which  typically  occur  when  the  sun 
heats  the  ground.  Buoyantly  stable  conditions  {C,>  0)  typically  occur  when  the  ground  cools  at 
night.  In  the  limit  -»  0 ,  the  T-functions  equal  0  and  the  profiles  assume  their  familiar 
logarithmic  forms  for  neutral  conditions. 

Now,  for  the  effective  sound  speed,  we  have 
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Using  the  thermodynamic  equations  Cp  =  Cv  +  R  and  y  =  Cp/Cv,  along  with  the  ideal  gas  law 
cl  =  yRT0 ,  we  find  Yd  =  gT0  (y  -  l)/c02 .  This  relationship  allows  the  second  term  on  the  right  of 

the  previous  equation  to  be  rewritten  as  g(y -l)(z -zr)/2cl .  Following  Wilson  [16],  it  is  also 
convenient  to  define  c ,  =  c0TJ2T0 .  We  then  have 
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Let  us  next  define  the  following  nondimensional  variables:  ceff  =  ce(f/c0 ,  c,  =  cjc0 ,  u,  =  ujc0 , 
8  =  8  /(c0 /) ,  z  =  zf/c0 ,  z,  =  zvg!cl ,  and  z0  =  z0g/c20 : 
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The  profile  cefr(z /  therefore  depends  on  the  nondimensional  parameters  z,.,  z0,g,  ct,  ut,  cos  /?. 
It  may  initially  seem  peculiar  to  normalize  zr  and  zo  with  cq  /  g  instead  of  co/f.  as  was  used  for 
other  lengths.  The  advantage  of  the  cq  /  g  normalization  is  that  zr  and  zo  are  fixed  values  in  many 
applications;  therefore,  zr  and  z0  are  constants  and  need  not  be  used  as  parameters  in  the  neural 
net. 


3.  Results 


The  development  of  these  nondimensional  equations  results  in  13  nondimensional  parameters 
that  need  to  be  specified.  These  are 

1 .  zs ,  the  source  height 

2.  zc ,  the  receiver  height 

3 .  xs ,  the  horizontal  separation  between  the  source  and  receiver 

4.  a  ,  the  normalized  static  flow  resistivity  for  the  ground 

5.  q  ,  the  tortuosity  to  porosity  ratio  for  the  ground 

6.  Na,  the  number  of  wavelengths  in  the  upper  absorbing  layer 

7.  Ji ,  the  attenuation  parameter  for  the  upper  absorbing  layer 

8.  zr ,  the  reference  height  at  which  the  temperature  is  measured 

9.  z0 ,  the  roughness  length 

1 0.  g  ,  the  normalized  gravitational  acceleration 

11.  ,  the  scale  for  turbulent  sound  speed  fluctuations 

12.  u, ,  the  scale  for  turbulent  wind  speed  fluctuations  (the  friction  velocity) 

13.  cos/?,  the  cosine  of  the  angle  between  the  propagation  and  wind  directions 
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Many  of  these  parameters  would  not  be  needed  for  the  training  of  a  neural  net.  If  the  calculations 
are  done  properly,  the  parameters  for  the  upper  absorbing  layer,  Nx  and  Ji ,  should  not  affect  the 
solution.  The  reference  height  zr  is  typically  a  constant  value.  For  a  given  ground  surface,  z0 , 
q ,  and  the  ratio  cr/g  are  constants.  57 g  is  a  constant  because  the  source  frequency  is  eliminated 
when  this  ratio  is  taken,  resulting  in  the  necessity  of  only  one  of  these  parameters  being 
specified.  Thus,  as  few  as  seven  parameters  could  be  used  to  train  the  neural  net.  These  include 
the  three  parameters  that  specify  the  source  and  receiver  locations  (zs,zc,  and  xs ),  the  three 

parameters  needed  to  calculate  the  normalized  effective  sound  speed  (ct ,  «, ,  and  cos/?),  and  the 
nondimensional  flow  resistivity  (a),  which  is  used  to  calculate  the  nondimensional  ground 
impedance. 

In  order  to  get  an  idea  of  the  relative  importance  of  each  of  these  seven  parameters  in  generating 
a  training  algorithm  for  an  artificial  neural  network,  it  would  be  necessary  to  determine  their 
sensitivity  on  the  sound  pressure  attenuation  at  a  specific  location.  This  was  performed  for  the 
seven  parameters,  with  the  results  indicated  in  Figures  1  and  2.  The  parameter  representing  each 
respective  curve  in  these  figures  was  varied  over  a  typical  physical  range,  and  the  resulting  sound 
pressure  attenuation  was  determined  from  the  nondimensional  model.  The  smallest  and  largest 
values  for  each  parameter  represent  0%  and  1 00%,  respectively,  in  the  figures.  These  values  and 
the  baseline  values  for  the  model  runs  in  which  the  respective  parameter  was  not  varied  are  given 
in  Table  1.  For  simulation  runs  in  which  the  receiver  height  and  horizontal  separation  were  not 
the  parameters  being  varied,  the  sound  pressure  attenuation  was  obtained  at  a  location  for  which 
the  value  of  these  parameters  were  1  m  and  500  m,  respectively. 

The  effects  of  varying  the  geometric  parameters  are  indicated  in  Figure  1  while  the  influence  of 
the  parameters  used  to  calculate  effective  sound  speed  and  ground  impedance  are  shown  in 
Figure  2.  Although  the  range  of  the  sound  pressure  attenuation  is  different  in  these  two  figures, 
their  scales  are  equivalent.  This  allows  a  direct  comparison  of  the  sensitivity  of  the  sound 
pressure  attenuation  to  each  of  the  seven  manipulated  parameters.  As  indicated  in  these  figures, 
the  attenuation  is  relatively  insensitive  to  the  source  and  receiver  heights  and  friction  velocity. 

For  the  source  and  receiver  heights,  this  is  attributable  to  the  small  range  over  which  these 
parameters  were  manipulated.  In  general,  varying  static  flow  resistivity  had  little  effect  on  sound 
attenuation.  For  small  values  of  this  parameter,  however,  there  was  a  more  substantial  change  in 
the  sound  pressure  attenuation.  This  could  be  partly  attributable  to  the  large  range  over  which 
this  parameter  was  varied,  which  was  almost  two  orders  of  magnitude.  As  the  wind  direction  was 
varied  from  0  to  180  degrees  (i.e.,  from  the  same  to  the  opposite  direction  of  sound  propagation), 
the  sound  attenuation  varied  by  approximately  1 0  dB.  The  ranges  over  which  the  normalized 
horizontal  separation  and  sound  speed  fluctuation  scale  were  varied  had  the  largest  effects  on 
attenuation,  with  both  of  them  producing  more  than  a  1 7-dB  reduction  in  the  sound  pressure.  The 
effect  of  the  sound  speed  fluctuation  scale  on  attenuation  was  unique  in  that  it  was  the  only 
parameter  that  did  not  exhibit  a  monotonically  increasing  or  decreasing  relationship. 

Because  the  sound  pressure  attenuation  values  are  normalized  to  pQc\ ,  it  is  probably  necessary 

to  put  the  results  in  proper  perspective.  Dividing  the  sound  pressure  by  the  product  of  p^c\ 

results  in  a  decrement  of  approximately  103  dB;  thus,  one  would  need  to  add  103  dB  to  the 
values  indicated  in  Figures  1  and  2  to  obtain  the  actual  sound  pressure  level  attenuations  in  the 
figures. 
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Figure  1.  Influence  of  normalized  source  height,  receiver 
height,  and  horizontal  separation  between  them 
on  the  normalized  sound  pressure  attenuation. 
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Figure  2.  Influence  of  the  angle  between  propagation 

and  wind  directions  (/?),  and  normalized  friction 
velocity,  scale  for  turbulent  sound  speed  fluctuations, 
and  static  flow  resistivity  on  the  normalized  sound 
pressure  attenuation. 
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Table  1.  Range  and  baseline  values  for  the  seven  normalized  parameters  (Actual  [dimensional] 
values  and  their  units  are  indicated  in  parentheses.) 


Parameter 

smallest  value 

largest  value 

baseline  value 

source  height,  zs 

0.0 

0.74 

0.148 

(m) 

(0.0) 

(5.0) 

(1.0) 

receiver  height, 

0.0 

0.74 

0.148 

(m) 

(0.0) 

(5.0) 

(1.0) 

horizontal  separation,  xs 

14.8 

133.4 

74.1 

(m) 

angle  between  propagation 

(100) 

(900) 

(500) 

and  wind  direction,  p 

0 

K 

0 

friction  velocity,  w* 

1.48  *10'4 

1.48  *10"3 

8.32*10  4 

(m/s) 

turbulent  sound-speed 

(0.05) 

(0.5) 

(0.28) 

scale,  c* 

3.5  *10'3 

0.296  *10~3 

-0.593*10“3 

(m/s) 

(-1.18) 

(0.1) 

(-0.2) 

static  flow  resistivity,  a 

200 

16000 

1605 

(mks  rayls/m) 

(25000) 

(2000000) 

(200000 

4.  Summary 


A  system  of  nondimensional  equations  for  sound  propagation  in  the  near-ground  atmosphere  was 
developed.  They  are  based  on  the  narrow  angle  parabolic  equation  for  sound  propagation  in  a 
moving  medium,  a  relaxation  model  for  the  sound-absorbing  lower  boundary  condition,  and 
Monin-Obukhov  similarity  for  the  atmospheric  profiles. 
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